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1 Introduction 

The primordial inflation is a well tested paradigm for the early universe [1], which is re- 
sponsible for creating the perturbations and the matter in the universe, for recent reviews 
see Refs. [2, 3]. Typically inflation can happen via slow rolling of a single scalar field on 
a smooth potential, which has unique predictions for the cosmic microwave background ra- 
diation [4]. The induced perturbations are generically random Gaussian fluctuations with 
almost scale invariant spectrum with a small tilt which indicates that inflation must come to 
an end in our patch of the universe. The fluctuations are known to be adiabatic perturbations 
which sources directly the curvature perturbations relevant for the structure formation. 

However, if there were two or more scalar fields rolling simultaneously on top of the 
potentials, there would be more sources for the curvature perturbations. It has already been 
known for a while that a relative perturbation from other fields, sometimes known as the 
entropy perturbations can also source the curvature perturbations along with the collective 
adiabatic perturbations during inflation [4-6]. A non- negligible coupling between the fields 
will typically curve the infiationary trajectory and the perturbations can be projected tangen- 
tial to the trajectory (adiabatic perturbations), and perpendicular to the trajectory (entropy 
perturbations) [6]. 

In order to track all the components of the perturbations, we base our formalism on 6N 
approximation (where A^ being the number of e-foldings) and separate universe approach [7- 
10]. The 6N formalism has been employed successfully for single and multi-field models of 
inflation [11]. Attempts to deal with the simplest forms of two-field potential — the separable 



potential by sum and product — have been carried out in Refs. [12-15] with the 5N forniahsm 
proposed in Refs. [7, 16]. 

In Ref. [12] when deahng with the separable potentials by product, the authors have 
suggested decomposing the perturbations into the adiabatic and entropy directions. Such 
a choice of decomposition has been deployed by many follow-up papers [13-15]. The final 
formulae for the relevant observables are usually complicated and indirect and usually the 
adiabatic limit is taken, see Refs. [13, 17, 18]. There have been attempts recently on predicting 
cosmological observables using moment transport equations, see Ref. [19-22]. Their results 
are still complicated enough and no analytical expressions are available for non-separable 
potentials. 

In this paper we decompose field perturbations into curvature and isocurvature per- 
turbations, which are defined differently from adiabatic and entropy modes. This allows us 
to present analytical solutions for separable and non-separable inflaton potentials. Our for- 
malism provides straightforward formulae for separable potentials, and the transparency in 
implementing the numerics for the non-separable inflaton potentials. Furthermore, it can be 
employed to study the bi-spectrum and tri-spectrum during inflation. 

It is well known that single field models of slow roll inflation do not generate large non- 
Gaussianity [23]. One requires a significant contribution from entropy perturbations arising 
from multi fields, and the conversion from entropy into curvature perturbations. The latter 
needs to be non-adiabatic in nature, for instance the end of inflation must happen via phase 
transition. 

In order to generate large non-Gaussianity, the requirements for a non-trivial boundary 
condition and conversion of entropy perturbations were first stressed in the context of pre- 
heating [24-26] . The preheating is a phenomenon of non-adiabatic particle creation, after the 
end of inflation when the inflaton is coherently oscillating, for a review see [27]. It was shown 
that during preheating one can generate significant non-Gaussianity from the conversion of 
entropy perturbations [26]. 

One can realize a very similar scenario during multi-field inflation, where the end of 
inflation happens not by the violation of slow roll conditions but by phase transition due 
to a waterfall field, as shown in Refs. [28], in the context of multi-brid inflation. In these 
models the end of inflation boundary and the inflaton trajectory intersects at an acute angle 
in the field space which converts the entropy perturbations into curvature perturbations at 
the end of inflation, and is also responsible for enhancing the observed non-Gaussianity. In 
this paper we will provide two examples: one with a separable potential to demonstrate how 
it works, and one with non-separable potential to show its capability, both of which acquire 
large non-Gaussianities. 

Since we base the formalism on 6N approximation, it still inherits systematic errors 
coming from 6N formalism. Such errors typically come from neglecting the interference be- 
tween different modes and taking field perturbations as uncorrelated Gaussian at the Hubble 
exit [29]. They may also come from the assumption that modes freeze once they go super- 
Hubble [30] . We will not discuss the significance of these errors as they are beyond the scope 
of this paper. 

2 Formalism 

Let us first decompose a generic perturbation into perturbations along the adiabatic and 
entropy directions as shown in Fig. 1. The adiabatic perturbation is converted into the 
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Figure 1. This is a schematic figure of trajectories and perturbation decompositions in the field space. 
The blue curve is a general trajectory, on which a general perturbation (Sff), Sip) at a point is shown by 
the black arrow. It can be decomposed to perturbations of fields Scj) and S^ shown by green arrows, or 
to adiabatic and entropy perturbations in blue arrows, or to curvature and isocurvature perturbations 
in red arrows. S(f) is taken as an example shown in purple that any perturbation including field 
perturbations can be decomposed to a curvature component "cur^^ and an isocurvature component 



classical curvature perturbation at the Hubble exit, while the entropy perturbation can still 
be converted into the curvature perturbation after that. For instance, in hybrid models there 
is a possibility of converting entropy into curvature perturbation at the end of inflation due 
to phase transition, see [28, 31, 32]. For a curved trajectory entropy perturbation can also be 
converted into curvature perturbation during inflation after the Hubble exit. For this reason 
we do not take the isocurvature direction at the Hubble exit the same as the entropy direction. 
Instead, we define the isocurvature direction as the one along which field perturbations do 
not generate any curvature perturbations, i.e. on a uniform- A^ hypersurface. We take the 
convention of A^ as the remaining e-foldings till a uniform energy density hypersurface after 
which the conversion from entropy to curvature perturbation is negligible. This gives 



dN 



-Hdt, 



(2.1) 



where H = dlna/dt is the Hubble expansion rate and a is the scale factor. For simplicity 
we will always assume the slow roll conditions. 

Therefore, in this paper we will decompose the field perturbations into curvature and 
isocurvature instead of adiabatic and entropy perturbations. The curvature direction is 
always along the trajectory of inflation, so it is actually the same as the adiabatic direction. 
Since their other respective components have different directions, the decomposition leads 
to different amplitudes for adiabatic and curvature perturbations, as demonstrated in Fig. 
1. This difference in amplitudes is exactly the amount of curvature perturbation converted 
from the entropy perturbation after the Hubble exit. 

We would also like to mention here that the decomposition into curvature and isocurva- 
ture perturbations is only initially known at the end of inflation hypersurface. The decompo- 



sitions during inflation should then be calculated through the method that is followed. Such 
calculations are necessary because all the information needed for cosmological observables is 
encoded in the decomposition, see section 2.4. 

2.1 Conventions 

Let us consider a general two-field slow roll inflation with the Lagrangian density 



^9^09^<A-^' 



£ = --9^09^(/.--a'^V'9^^-y((/>,V), (2.2) 



where (j) and ip are the two slow roll scalar flelds whose perturbation distributions are con- 
sidered independently Gaussian, and 1^(0, V') is their potential. For convenience, we use a 
slightly different convention for the slow roll parameters 

''' = 47^' ^^-^^ 



W^ = -^^^^ (2-5) 

e^^el + el- (2.6) 

Here Mp = 1/VG is the the Planck mass, and the subscripts </> and V' after ",'' denote 
derivatives with respect to the fields. For simplicity, unless expressions for cj) and ip are 
asymmetric, we only give the equations for (/> here and onwards. 

As long as the quantum effect does not overwhelm the classical slow roll, we can define 
C{(j), ip) as the variable to track the trajectory whose exact expression can be found in Eq. 
(A. 6). In fact it can be chosen arbitrarily as long as it is a unique function of the trajec- 
tory, at least locally. It should also be defined as such that the infinitesimal isocurvature 
perturbations have nonzero dC{(j), ip) and vanishing 5N, while the curvature perturbations 
have non- vanishing 6N and vanishing 5C{(p,^), which indicate dN/dC = dC/dN = 0. As a 
result of this the curvature and isocurvature perturbations never transfer into each other. 

Therefore we have a two-dimensional phase space during infiation, represented by pa- 
rameters either (</>, V') or {N,C). Either parameter set can thus be expressed as a unique 
function of the other set, such as N(4>,ij;) and 4>{N,C). For a specific trajectory, C is a 
constant, so we may omit it sometimes and write (f){N) only. For a general partial derivative, 
we will use the variable explicitly in the subscript after "," ^■ 

2.2 Background and first order perturbations 

According to SN formalism, see Refs. [7, 9], on large scales calculating the power spectrum 
of curvature perturbation is equivalent to calculating the power spectrum of 6N in separate 
universes on the initially flat hypersurface at the Hubble exit. The uncorrelated Gaussian 
perturbations (50,5^) of (</>, V')) then generate the perturbations in the e- foldings: 

6N = N^^6(P + N^6^ + - (iV,0^(5</.2 + N^^St/;^'^ + N^^^Scf^di; + higher orders. (2.7) 



^Primes will denote the partial derivatives w.r.t TV while keeping C constant but regarding 4> and tp as 
functions of A'^, and dots will denote taking partial derivatives w.r.t C keeping A'' constant in the same way. 



The power spectra of 5(j) and 5ip take the same value: 

'Pscf} = 'Psil, = -7^- (2.8) 

Now consider the period when the slow roll conditions hold for both the fields. The 
equation of motion for (j) then becomes: 

_|+3i/^ + F, = 0. (2.9) 

After applying the slow roll approximations and Eq. (2.1), it reduces to 

d<^=!f^diV. (2.10) 



In a separate universe approach where super-Hubble perturbations are smoothened to a 
locally homogeneous perturbation {5(j),6'tp), the first order perturbation from Eq. (2.10) can 
be written as: 

d6^ = '-f^d6N + ^6e^dN. (2.11) 

2V7r 2-y/7r 

Since our aim is to calculate N^^p, we take Sep and 5^l^ to be infinitesimal. As shown in Fig. 1 
in purple, we decompose 6(j) in two directions: 

• Curvature direction: along the trajectory defined as 6(l)cur, and 

• Isocurvature direction: along the uniform- A^ hypersurface defined as 5(f)iso- 
According to this definition, their separate equations from Eq. (2.11) are 

d5(l},ur = ^e^dSN, (2.12) 

d5(t)iso = ir^5eA,dN 



2e^)6(j)iso + {'n<t>i> - 2e^e^)6i)iso)dN. (2.13) 

Note that these two equations have different meaning. Eq. (2.13) tells us how Scpiso evolves as 
the universe is evolving, while Eq. (2.12) tells us how many SN would be generated (linearly) 
by the perturbations 5<j)cur- So we can remove the d's in Eq. (2.12) and rewrite it as 

^^ - ^"-.e^SN. (2.14) 



2^AF 



Similar expressions can be derived for the if) field. The Eq. (2.13) together with the ip version 
becomes a set of differential equations for 6(j)iso and Sipiso, provided with the background 
trajectory, C. These differential equations have the set of solutions given by: 

'''' 'dC = U{N,C)6C, (2.15) 



dC 



where /^ is some function of A^ and C, and 5C is the deviation from the original trajectory 
which stays constant as the separate universe evolves. 



Since the curvature and isocurvature directions are not necessarily perpendicular, solv- 
ing N^fj, becomes more complicated. First of all, it is defined as 



ON 5N 

N,^=-—=—- 



(2.16) 

5ip=0 



+ 5<P.so= ''^i^'f* M^6N, (2.18) 



Constructing the 6ip = Stpcur + ^''Piso = direction from the V' version of Eq. (2.14) and Eq. 
(2.15) yields: 

5C = -^^^5N. (2.17) 

2v^/«/' 

This is followed by 



which provides 

Since e^ and /^ are both functions of N for a given trajectory, we solve N^^ by taking the 
derivate w.r.t A^ on both sides of Eq. (2.19), while (j) and V' are regarded as functions of A^. 
After some calculation and substituting Eq. (2.13), we reduce it to 

N[^ = U4,N^^ + v^,, (2.20) 

where the prime denotes taking derivative w.r.t N while fixing C, and 

U(i, = — jy^v^ - '^<^<A! (2-21) 



The above differential equation has an exact analytic solution for N^^ as a function of N (for 
a non- vanishing v^) 

N^{N) = iV,^(iVo) + e^^o"^(^)'i^ r v^iN)e-^^o''^^^^''''dN. (2.23) 

Jno 

Here N^ff){No) is the value of N^^p at N = Nq, which serves as the boundary condition ^. 

2.3 Second order perturbations 

In certain cases, such as in the case of separable potentials [12, 13], the integrals in Eq. (2.23) 
can be worked out analytically, so A^^^ becomes a simple function of (j) and ip. For such cases 
the second order derivatives A^,00, A^.^i/) and N^^^ can be derived by simply differentiating 
N^fj, and -/V,^, while keeping aware of nonzero Nq^^ for general end of inflation and/or post- 
inflationary conditions, such as in the case of reheating or preheating. In general, such 
integrals cannot be solved analytically. To calculate the second order derivatives for any 
potential, here we provide a general framework. 



^Although in Eq. (2.23) the integral of u^ niay give readers the iUusion that it can be worked out analyticaUy 
for any potential, it actually cannot because u^ is a function of both (f} and tp, neither of which remains constant. 
Therefore one cannot simply change the integral variable to ({> or tp only. 



By definition, we decompose the second derivative along (j) direction to the curvature 
and the isocurvature directions 

N^^^ = ^ = N[^N^^ + N^C,^. (2.24) 

Since we already know the first term on the r.h.s, we only need to calculate the second term. 
To solve N^^ we can take dots on Eq. (2.20) and utilize the interchangability between primes 
and dots. We then obtain the equation for A^,^, 

N'^ = u^N^^ + (u^^^N^^ + v</,,0)/0 + {u^^^N^^ + v^^^)f^. (2.25) 

Note that dot denotes derivative w.r.t C . This equation holds the same form as that of Eq. 
(2.20), and the solution should have the same integrated form as that of Eq. (2.23). In Eq. 
(2.24) and Eq. (2.25), C,^ and /^ can be solved by simple integration through 

C[^ = u^C^^, (2.26) 

4 = {r}<t><t> - 2e^ - §^(W^ - 2606^)) U- (2-27) 

We can also obtain N^^ by differentiating Eq. (2.23) w.r.t C. It is expected to have 
exactly the same solution as that of in Eq. (2.25), as long we keep in mind that Nq may 
also be a function of C (if the boundary is not on a uniform- A^ hypersurface) , so Nq / 0. 
Similarly, we can also calculate 

N^^^ = ^ = N[^N^^ + iV,^C^. (2.28) 

Although we are capable of presenting the full integrated form for A^.,/,,/,, we are not 
expressing this in the current paper for obvious reasons — the expressions are long and not 
practical for general non-separable potentials and we are not using the full integrated form in 
this paper. However, they do show that there exist analytical solutions for two-field slow roll 
inflation even for non-separable potentials albeit in the integrated form. In practice when 
calculating models with non-integrable N^^f, (in Eq. (2.23)), differential equations Eq. (2.20), 
Eq. (2.25), Eq. (2.26), and Eq. (2.27) are actually far more useful — we can combine them 
with the background equations of motion and obtain all of them in a single run by numerically 
solving the differential equations backwards from the end of inflation to the Hubble exit with 
an appropriate boundary condition. 

Regarding the boundary conditions, the derivation for a general end of inflation con- 
dition is given in the appendix section A. For separable potentials by sum or product, the 
integrals can be worked out analytically, and therefore the results can be simplified and 
extended. Such calculations are shown in the appendix section B. 

2.4 Cosmological observables and non-Gaussianity 

We now have all the parameters required for deriving the cosmological observables. The 
power spectrum of curvature perturbations is equal to the power spectrum of 5N on large 
scales, which means 

P^ = PsN = N^Ps^, + N%Ps^ = ^{N^ + N%). (2.29) 



The scalar spectral index, Ug, is defined as 



n. 



dlnp2 
dA^ 



-2e' 



2{N'^^N^ + N'^^N^) 



Nl + N 



(2.30) 



V> 



whose running can be determined by simple differentiation while using Eq. (2.20). The 
amplitude of bi-spectrum is characterized by the parameter /nl, see [16, 29]: 



/nl 
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(2.31) 



For tri-spectrum, we can only calculate up to second order in this formalism, which gives one 
of the two tri-spectrum parameters, see [33] 
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(2.32) 



We define the non-adiabaticity parameter a for field perturbations at the end of inflation, 

(2.33) 



a 



1 — a 






where P? and Pr indicate the power spectra of field perturbations in entropy and curvature 
directions respectively, and the subscript e denotes the end of inflation. The correlation 
parameter is defined from the cross-correlation power spectrum P??, as 
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P~sc 
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(2.34) 



We calculate the field perturbations along the entropy and the curvature directions at the 
end of inflation from those at the Hubble exit through 5N and 6C, which remain constant 
during inflation. This yields 
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(2.35) 
(2.36) 



where ± should be chosen as the sign of (e,^/^ — €^f^)\e- The definition of a and /3 are in the 
field space, unlike a and /3 defined as power spectrum ratios of energy density perturbations 
in Ref . [34] . This is because one cannot split the total energy density to field energy densities 
if they are coupled. However, for separable potentials by sum in the form Eq. (B.l), fields 
are uncoupled so we can calculate the power spectra of non-adiabaticity perturbation and its 



correlation with curvature perturbation in energy density. See the derivation in the appendix 
section B.l. 

We also define the power spectrum ratio between adiabatic and curvature perturbations: 

where 7 — )• means the whole curvature perturbation comes from the conversion of the 
entropy perturbation and vice versa. Strictly speaking the approximation of taking the 
adiabatic perturbation as the curvature perturbation (i.e. the so-called adiabatic limit) is 
only allowed when 7 — )• 1 is satisfied both at the time of Hubble exit and onwards. 

3 Two-field infiationary models 

We will demonstrate our formalism through a simple model of two-field inflation in which 
two scalar fields cj) and ip prevent the waterfall field x from falling ^. The potential is given 
by: 

y = yo(i-^] +{gy+glijV + vi4>,i^). (3.1) 

Infiation is dominated by the constant potential energy density, Vq, with x sitting at VEV 
during infiation. Both (p and ip are slow rolling initially ^. Inflation ends once the effective 
mass of X becomes negative, which gives the end of inflation condition as 

gl(t>l + 9l^I = 2Vo/v^ (3.2) 

where the subscript e means the end of inflation. 

The end of inflation happens by waterfall, so the trajectory of inflation and the end 
of inflation boundary can form a sharp angle, which allows a significant transfer from the 
entropy to curvature perturbations at the end of inflation. Large non-Gaussianity can also be 
achieved through this mechanism when the couplings g^j) and 5^ differ significantly. This also 
requires that the entropy perturbation to remain non- vanishing during the last 60 e-foldings 
of exponential expansion, which means one field has to be lighter than the other. This is 
shown in Fig. 2. 

Since during inflation the effective potential is just V = Vq + V{<j),^l'), whether the 
potential is separable depends on V. In the following, we will first demonstrate a separable 
case, giving both analytical and numerical results for comparison, and then provide a non- 
separable potential model. 

3.1 Separable potentials 

For the separable potential case, we take the double quadratic potential of the form: 

V{^,%b)=ml<P^+ml^\ (3.3) 

^Similar studies were performed in Refs. [28, 32, 35-38]. Note that in this paper (f> and ip are not charged 
under the Standard Model gauge group or its minimal extensions [39-41]. Here we treat them as singlets, and 
we also assume that 0, ip and x all decay into the Standard Model radiation without making any justification. 
In this respect this model serves at best as a simple toy model which can only mimic the success of CMB 
perturbations. 

*We do not study the initial condition problem for this model, we expect that the initial condition problem 
for this model will be similar to that of any hybrid model of inflation, see [42]. 



Isocurvature Entropy 




Figure 2. A schematic figure shows how the entropy perturbations are converted into the curvature 
(and isocurvature) perturbations at the end of inflation, by assuming that no significant conversion 
happens after inflation. For Sip > 0, the trajectory is shifted upwards slightly, leading to a longer 
inflation and therefore SN > 0. A large non-Gaussianity may be generated if the boundary is very 
curved at its intersection with the trajectory. 



We then use the expression of A^^,^ for separable potential case by substituting Eq. (3.3) into 
Eq. (B.2). We also take the approximation that ip is much lighter than (p, so it moves very 
slowly during the last 60 e-foldings of inflation. The dominating terms, at first order are: 
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AttVo 
M^miiP 



(3.4) 
(3.5) 



where 9 characterizes the end point of inflation on the elliptical boundary, and it is defined 
as: 



tan^ 



9i^1pe 



(3.6) 



For the case rn'^/m'^ ^ tan'^ 9 which we are interested in, the differential relations hold: 



89 _ rn^ tan ( 
dcp m?i (j) 



89 _ tan 6* 

dip ip 



(3.7) 



The second order derivatives can then be calculated through differentiating Eq. (3.4) and Eq. 
(3.5), giving the dominant term 
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Figure 3. Red and blue curves are analytical and numerical results respectively. They are too 
close to be distinguished from each other, showing the analytical approximations as compared to the 
numerical results for the physical parameters Eq. (3.13). 

The cosmological observables can then be derived from substituting Eq. (3.4), Eq. (3.5) and 
Eq. (3.8) into Eq. (2.29), Eq. (2.31) and Eq. (2.37), giving 
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m\ + rnL tan^ 9)"^ 



mh + ml tan^ i 



+ 7T1-1 tan^(0)(/)2/^2 + m\ilj'^ / (fP' 
The slow roll approximation gives the evolution of fields 



=,e '^''^0 
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(3.9) 

(3.10) 
(3.11) 

(3.12) 



Since m^ <^ m^, during inflation ip remains almost constant and (p drops gradually, 
while all other parameters remain constant. Therefore as N gradually decreases, we obtain 
an increasing trend of P^ (which remains within 2a of the current bound on the amplitude 
of the power spectrum [1]), and a decreasing /nl, while a slowly increasing n^ — 1, and 7 
depending on the parameter choice (which is increasing here). 

Our formalism allows us to numerically solve the differential equations Eq. (2.20), Eq. 
(2.25), Eq. (2.26), and Eq. (2.27). Unlike the finite-difference method applied in Ref. [13], this 
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Figure 4. Numerical results of bi- and tri-spectrum parameters with fixed model parameters in Eq. 
(3.13), but varying initial conditions of (f) and ip. See Eq. (3.6) and Eq. (3.11). 

method is not plagued by the dilemma of balancing between error sources ^, so it typically 
has much more precise results at the second order and the errors are easily controllable. 

For this model, numerical results are shown for the parameters in Eq. (3.1) and Eq. 
(3.3): 



m^ = 3.03 X 10^° GeV, 
m^ = 9.57 X 10^ GeV, 



9^} 



1.13 X lO"^ 

1.79 X 10^1, 



V = 



= (1.09 X 10^^ GeV)^, 
1.22 X 10^6 GeV. 



(3.13) 



For these parameters if the fields start from -0 <C </>, we can obtain a large non-Gaussianity. 
As an example, here we consider cj) = 1.95 x 10^^ GeV and ip = 3.91 x 10^^ GeV, which result 
in the total scalar power spectrum, scalar spectral tilt, non-Gaussianity and the adiabatic- 
curvature ratio as functions of remaining e-foldings shown and compared with the analytical 
solutions in Fig. 3. In this case, we also obtain the largest tnl ~ 20000. 

It is worth mentioning here that if we change the initial condition while fixing all the 
parameters, we are able to get /nl and tnl as functions of 7 at the pivot scale, corresponding 
to A; = 0.002 Mpc~ , which are shown in Fig. 4. We can see that under these parameters 
an interesting result is that we obtain the largest /nl when 7 = 0.5, i.e. when adiabatic 
and entropy perturbations give equal contributions to the power spectrum of the curvature 
perturbations. 

This effect is easily understood qualitatively. A dominant contribution to /nl comes 
from the entropy perturbations so we expect /nl to grow as 7 decreases from 1. As long as 
changing 7 does not significantly affect the value of (j) at the Hubble exit of the relevant scale, 
the relation N^ -\- N^ oc 7"^ holds which contributes a 7^ coefficient to /nl- Also, one can 
find the same conclusion by solving as a function of 7 and substituting it into /nl > which 
then precisely leads to /nl oc 7(1 — t)- For tnl, a similar relation holds tnl o<^ 7^(1 — t)- 

^In the finite difference method errors come from the inaccuracies from higher order contributions and 
imprecisions from small relative differences, whose strengths are inversely correlated. If the pivot scale exits 
the Hubble patch during infiation leaving /nl ~ 100 and (j)N^4, — 10^A''pivot, where A^'pivot ~ 50 has 10 
significant figures from solving the equations numerically, then a rough estimation of the error for uncorrected 
finite difference method gives the number of significant figures of A^,</></, at most 1.75 — l2:|/4. 
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Figure 5. Numerical results for a non-separable potential is shown above. For Fig. 5(a) and Fig. 
5(b), we use the fixed initial conditions. For Fig. 5(c) and Fig. 5(d), we allow the initial conditions 
of (j) and ^ to vary while fixing the e-folding we are interested in. The physical parameters are given 
by Eq. (3.15). 



3.2 Non-separable potentials 

As a non-separable example, we consider a simple Logarithmic potential 



V{(t),ij) = Vi\n 



M2 
p 



(3.14) 



where Vi is the energy scale of the two fields, and A is the parameter which characterizes the 
mass ratio between them. 

For the parameters shown below, we obtain the right amplitude and the tilt in the power 
spectrum, but also large non-Gaussianity 



(/> = 1.59 X 10^6 GeV, 
V^ = 7.18 X 10^° GeV, 
V = 1.22 X IQi'^ GeV, 



04, = 3.16 X 10-5, 
5^ = 1.30 X 10-2, 
A2 = 0.1. 



Vq = (1.83 X IQi^ GeV)^ 

Vi = (2.58 X 10^2 GevK (3.15) 



The numerical solutions of Eq. (2.20), Eq. (2.25), Eq. (2.26), and Eq. (2.27) are shown in Fig. 
5, with P(^ ~ 2.4 X IQ-^ and n^ — 1 w — 1.3 x 10-^, which are quite similar to Fig. 3 and Fig. 
4. This is reasonable as they rely on the same end-of-inflation mechanism to generate large 
non-Gaussianity. From this example, it is evident that the numerical method also works well 
with non-separable potentials. 
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We note that the /nl's and tnl's have similar shapes in Fig. 4 and Fig. 5. The large 
non-Gaussianity arises due to the end of inflation boundary condition rather than the insep- 
arability of the potential. Although here the inseparability is not playing an important role 
and we only use it to demonstrate the capability of the formalism, non-separable potentials 
still provide distinct features in many cases, such as multi-stream inflation[43] and the study 
of the effects of couplings for N-flation[44, 45]. In such models this formalism will be helpful 
in the prediction of cosmological observables. 

4 Comparison with related works 

At first, Yokoyama et al. proposed to calculate the evolution of N^^ and A^,^;/ backwards 
from the end of inflation to the Hubble exit in Refs. [46, 47]. They used the transfer matrix 
method, and the analytical representations are mostly formal. In particular, their derivation 
of N^^i, has two levels of embedded integrals which are unfavored by numerical calculations. 
Recently Mulryne et al. managed to evolve the distribution of fleld perturbations after their 
Hubble exit, which is called moment transport equation method [19, 20]. These methods are 
physically equivalent but have different advantages — in a single run the backward formalism 
is capable of calculating the amplitudes of all modes at a specific time, which is the case for 
CMB, while the moment transport equation method evolves a specific mode all through 
inflation, deriving its amplitude at any time. 

This paper has the same idea with those by Yokoyama et al, but attempts to understand 
the perturbations in a geometrical approach. Therefore the derivations are more visual 
and straightforward than Refs. [46, 47]. In addition, two improvements have been made 
in this paper. As a result of the geometrical approach and the additional parameter C, 
we are able to address isocurvature perturbations within the framework of 5N formalism, 
which is otherwise considered difficult as mentioned in Ref. [22]. We also use a generic 
end of infiation condition, which allows post-infiation and/or end-of-infiation mechanisms 
to generate significant curvature perturbations. For two-field inflation, we also find the 
equations for derivatives of N (e.g. Eq. (2.23) and Eq. (2.25)) can be separated, providing a 
higher efficiency in cases where we know in advance whose perturbation will dominate. 

We would also like to briefly mention here one specific advantage of the backward 
formalism (e.g. this paper) over moment transport equation method. As shown in section 
2.2 and section 2.3, derivatives of N and C are independent and can be solved separately. 
This means we can pick out only the one(s) we are interested in and save time by discarding 
the rest, which however the moment transport equation method is unable to achieve. The 
lower time complexity will become obvious when we have many fields such as in N-flation.^ 

We have compared the analytical differential equations of N's derivatives with those in 
Refs. [19, 20, 22, 46, 47], and they match in every detail. Verifications have also been made 
through comparisons between papers about separable potentials and our reduced result in 
section B. 



^For n-field calculation, Refs. [46, 47] introduced O^ to reduce the number of integrals to 0{n). Although 
the n-field extension of this formalism has 0{n^) integrals, the total time complexity remains the same with 
that of [46, 47]. This is because each integral in [46, 47] needs 0{n^) time to calculate every sum in each 
integral while (the extension of) this formalism only needs 0{n). Therefore both methods have a total 0{n^) 
time complexity and Refs. [46, 47] don't have any advantage over this formalism in this sense. 
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5 Summary 

We have proposed a simple framework based on 6N formalism which deals with two-field slow 
roll inflation models. The formalism gives analytical integrated solutions up to second order 
perturbations for general non-separable potentials. In the case of separable potentials, they 
are exactly integrable. This formalism can be easily adapted for numerical solutions to non- 
separable potentials. In the examples we demonstrated, our results match the expectations of 
[26, 28, 32] that in order to obtain large /nl and tnl, one must generate entropy perturbations 
during inflation which then convert to curvature perturbations non-adiabatically. 

The major difference with previous studies is that we decompose the perturbations into 
curvature and isocurvature rather than adiabatic and entropy. Therefore it also enables the 
calculation of isocurvature perturbations while basing on SN formalism. Although this paper 
only discusses two-field slow roll inflation, our formalism can also be extended to arbitrary 
number of fields with non-canonical kinetic terms, and also fast roll fields. 
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A Boundary conditions 

Let us suppose that the two-field slow roll inflation stage ends at some general classical 
boundary, given by: cr((/), ip) = 0, which could be either due to violation of slow roll conditions, 
such as (T = e^ — 1, or via some sudden phase transitions. Since the boundary need not 
necessarily be on a uniform- A'^(/isocurvature) hypersurface, we will assume on the boundary 
the remaining e- foldings Nq[(J)^iP) is well-known, which has dNo/dC ^ 0. 

From Fig. 6, we can work out the first and second order boundary conditions through 
vectors. Pick a point (0, ip) on the boundary and set C = here, the unit vector along the 
trajectory is then given by 

n.-^^ML = f^,^V (A.1) 

Along the boundary, we have 5a = a^,f,64) + a^^Sip = 0, so its unit vector is 

n,^^^£l^£^. (A.2) 

""% + ''I 

To find a vector on the uniform- A^ hypersurface, we construct an infinitesimal general vector 

X = xihi + XbUh. (A.3) 

This displacement in the phase space then generates 



6N = {Nq^^, A^o,^) • UbXb + '^xi, (A.4) 
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Trajectory 




Figure 6. A schematic figure which illustrates how C is defined through uniform- A^ hypersurfaces 
at the end of infiation. The blue and red curves are inflationary trajectory and the end of inflation 
boundary, respectively. The arrows in blue and red are unit vectors n; and n;,, from which the 
green arrow hat, the vector in the uniform- A^(/isocurvature) direction, can be constructed. The green 
dashed lines show the uniform- A^ hypersurfaces. 

in which the first term on the r.h.s is the contribution outside(/after crossing) the boundary, 
and the second term comes within the boundary (during two-field slow roll inflation). 

By vanishing 5N in Eq. (A. 4), we are able to find the vector on the uniform- A^ hyper- 
surface 



rvN = rife 



2j7:{al + a%) 



-n/. 



For a general displacement at the end of inflation, x = {5(j),5ip), we define C (or 5C) by 
decomposing it as 

y. = 5C\iN + xm\]si=NQ- (A.6) 

Then obviously the condition dN/dC = dC /dN = is automatically satisfied. From Eq. 
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(A. 6), we find the following boundary conditions by vanishing xi, Sip, or 6(1) 

^\/<^% + ^) 

CANo) = -^^ -, (A.IO) 

^ /^ ^ 2V7ra^/Mp - (A^p.^ay, - A^o,v,o-,0)e^ ^^ 

C,^{No) = ^4 -■ (^-12) 

Boundary conditions for the second order perturbations can be derived from 

iV.^A^o) = U{No)N^^{No) + UiNo)N^^{No). (A.13) 

B Separable potentials 

B.l Separable by sum 

Using the potential, 

V{ct),^P) = U{(|)) + W{^l)), (B.l) 

we will find the integral in Eq. (2.23) can now be worked out analytically 

N^ = ^^N^No) + ^{U- UiNo)) . (B.2) 

Eq. (2.26) also gives 

C^=^^^C4No). (B.3) 

In the above equations, we have used primes on separate potentials as partial derivatives 
w.r.t the field, Nq as A^ at the boundary, and U{X) as the value of C/ at A^ = A on a specific 
trajectory. After simple manipulation, we find Eq. (B.2) is exactly the same with the results 
obtained in Ref. [13]. For the second order derivatives, we can just take derivatives from Eq. 
(B.2), while paying extra attention on A'o,^, which may be non- vanishing. 

With the separable potential by sum Eq. (B.l), if the energy densities of (p ^-nd ^|J 
dominate over other components, we are able to calculate the non-adiabaticity in energy 
density perturbations. Note that for the separable potentials, mentioned in section 3, the 
total energy density is dominated by the vacuum energy density of the waterfall field Xi so 
the following discussion is not applicable to those cases. 
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For a separable potential, at the end of inflation we can define the non-adiabaticity 
perturbation in energy density, as 



S 



P<t> P4> 



(B.4) 



in which p^ is the energy density of (j) and 5p0 denotes its perturbation. Here the subscript 
e also means the end of inflation, i.e. taking the expressions at A^ = Nq. Under the slow roll 
approximations, we find 



yijU ^ /^ 



{C^^5<p + C^^5il:). 



(B.5) 



Remembering that C, = N^^5(f)+N^^5'il^, we obtain the power spectra of S and cross-correlation 
between S and (", as 



Ps = is') = (^/, - 




(C,9i + C,v)4^2' 


(B.6) 


Psc^{SO= (jtU- 


W' \ 


^iC^N^ + C,^N^)^. 


(B.7) 


Their spectral tilts are given by: 






T dlnP5 „ 2 
"^ ^ = dN = '' - 




(B.8) 


^^^PS( „ 2 

"^^ ^= dN = '' 


c<pn;^ 




(B.9) 



Following Ref. [34], we define the non-adiabaticity parameter a and the cross-correlation 
parameter /3 as 



a _ Ps^ 
l-a~ Pc_ 



fU' W 



/3 = 



P 



sc 



VP^C 






(B.IO) 
(B.ll) 



When comparing with the non-adiabaticity in the field perturbations, Eq. (2.35) and Eq. 
(2.36), we find that although the expressions of a and a have their boundary parts different, 
/3 and /5 only differ by a sign. 

B.2 Separable by product 

For separable potentials by product, we use the generic form 

V{(j), il^) = M^e^(<^)+^('^^ (B.12) 

We can easily obtain: 

The second order derivatives can be derived in the same way from taking derivatives from 
Eq. (B.13). 



18 



References 

[1] WMAP Collaboration Collaboration, E. Komatsu et. ai, Seven-Year Wilkinson Microwave 
Anisotropy Probe (WMAP) Observations: Cosmological Interpretation^ Astrophys.J.Suppl. 192 
(2011) 18, [arXiv: 1001.4538]. 

[2] A. Mazunidar and J. Rocher, Particle physics models of inflation and curvaton scenarios, 
Phys.Rept. 497 (2011) 85-215, [arXiv: 1001 .0993]. 

[3] A. Mazumdar, The origin of dark matter, matter- anti-matter asymmetry, and inflation, 
arXiv: 1106.5408. 

[4] V. F. Mukhanov, H. Feldman, and R. H. Brandenberger, Theory of cosmological perturbations. 
Part 1. Classical perturbations. Part 2. Quantum theory of perturbations. Part 3. Extensions, 
Phys.Rept. 215 (1992) 203-333. 

[5] H. Kodama and M. Sasaki, Cosmological Perturbation Theory, Prog.Theor.Phys.Suppl. 78 
(1984) 1-166. 

[6] C. Gordon, D. Wands, B. A. Bassctt, and R. Maartens, Adiabatic and entropy perturbations 
from inflation, Phys.Rev. D63 (2001) 023506, [astro-ph/0009131]. 

[7] M. Sasaki and E. D. Stewart, A General analytic formula for the spectral index of the density 
perturbations produced during inflation, Prog.Theor.Phys. 95 (1996) 71-78, 
[astro-ph/9507001]. 

[8] M. Sasaki and T. Tanaka, Superhorizon scale dynamics of multiscalar inflation, 
Prog.Theor.Phys. 99 (1998) 763-782, [gr-qc/9801017]. 

[9] D. Wands, K. A. Malik, D. H. Lyth, and A. R. Liddle, A New approach to the evolution of 
cosmological perturbations on large scales, Phys.Rev. D62 (2000) 043527, [astro-ph/0003278]. 

[10] M. Sasaki, Nonlinear curvature perturbations in an exactly soluble model of multi- component 
slow-roll inflation, Class. Quant. Crav. 24 (2007) 2433-2438, [astro-ph/0702182]. 

[11] A. R. Liddle, A. Mazumdar, and F. E. Schunck, Assisted inflation, Phys.Rev. D58 (1998) 
061301, [astro-ph/9804177]. 

[12] J. Garcia-Bellido and D. Wands, Metric perturbations in two field inflation, Phys.Rev. D53 
(1996) 5437-5445, [astro-pli/9511029]. 

[13] F. Vernizzi and D. Wands, Non-gaussianities in two-field inflation, JCAP 0605 (2006) 019, 
[astro-ph/0603799] . 

[14] K.-Y. Choi, L. M. Hall, and C. van de Bruck, Spectral Running and Non-Gaussianity from 
Slow-Roll Inflation m Generalised Two-Field Models, JCAP 0702 (2007) 029, 
[astro-ph/0701247]. 

[15] Q.-G. Huang, Negative spectral index of f^L in the axion-type curvaton model, JCAP 1011 
(2010) 026, [arXiv: 1008.2641]. 

[16] D. H. Lyth and Y. Rodriguez, The Inflationary prediction for primordial non-Gaussianity, 
Phys.Rev. Lett. 95 (2005) 121302, [astro-pli/0504045]. 

[17] J. EUiston, D. J. Mulryne, D. Seery, and R. Tavakol, Evolution of fNL to the adiabatic limit, 
JCAP 1111 (2011) 005, [arXiv: 1106.2153]. 

[18] J. Elliston, D. Mulryne, D. Seery, and R. Tavakol, Evolution of non-Gaussianity in multi-scalar 
field models, Int.J.Mod.Phys. A26 (2011) 3821-3832, [arXiv: 1107.2270]. 

[19] D. J. Mulryne, D. Seery, and D. Wesley, Moment transport equations for the primordial 
curvature perturbation, JCAP 1104 (2011) 030, [arXiv: 1008.3159]. 

[20] D. J. Mulryne, D. Seery, and D. Wesley, Moment transport equations for non-Gaussianity, 



19 



JCAP 1001 (2010) 024, [arXiv: 0909. 2256]. 

[21] J. Frazer and A. R. Liddle, Multi-field inflation with random potentials: field dimension, 
feature scale and non-Gaussianity, arXiv: 1111 .6646. 

[22] D. Seery, D. J. Mulryne, J. Frazer, and R. H. Ribeiro, Inflationary perturbation theory is 
geometrical optics in phase space, arXiv: 1203 .2635. 

[23] J. M. Maldacena, Non- Gaussian features of primordial fluctuations in single field inflationary 
models, JHEP 0305 (2003) 013, [astro-pli/0210603]. 

[24] K. Enqvist, A. Jokinen, A. Mazumdar, T. Multamaki, and A. Vaihkonen, Non-Gaussianity 
from preheating, Phys. Rev. Lett. 94 (2005) 161301, [astro-ph/0411394]. 

[25] K. Enqvist, A. Jokinen, A. Mazumdar, T. Multamaki, and A. Vaihkonen, Non-Gaussianity 
from instant and tachyonic preheating, JGAP 0503 (2005) 010, [hep-ph/0501076]. 

[26] A. Jokinen and A. Mazumdar, Very large primordial non-gaussianity from multi-field: 
application to massless preheating, JGAP 0604 (2006) 003, [astro-ph/0512368]. 

[27] R. AUahverdi, R. Brandenberger, F.-Y. Cyr-Racine, and A. Mazumdar, Reheating in 

Inflationary Gosmology: Theory and Applications, Ann. Rev. Nucl. Part. Sci. 60 (2010) 27-51, 
[arXiv: 1001.2600]. 

[28] M. Sasaki, Multi-brid inflation and non-Gaussianity, Prog. Theor. Phys. 120 (2008) 159-174, 
[arXiv: 0805. 0974]. 

[29] D. Seery and J. E. Lidsey, Primordial non-Gaussianities from multiple-field inflation, JGAP 
0509 (2005) Oil, [astro-ph/0506056]. 

[30] E. Nalson, A. J. Christopherson, I. Huston, and K. A. Malik, Quantifying the behaviour of 
curvature perturbations during inflation, arXiv: 1111 .6940. 

[31] A. D. Linde, Axions in inflationary cosmology, Phys. Lett. B259 (1991) 38-47. 

[32] A. Naruko and M. Sasaki, Large non-Gaussianity from multi-brid inflation, Prog. Theor. Phys. 
121 (2009) 193-210, [arXiv:0807.0180]. 

[33] C. T. Byrnes, M. Sasaki, and D. Wands, The primordial trispectrum from inflation, Phys. Rev. 
D74 (2006) 123519, [astro-ph/0611075]. 

[34] WMAP Collaboration Collaboration, E. Komatsu et. al., Five-Year Wilkinson Microwave 
Anisotropy Probe (WMAP) Observations: Cosmological Interpretation, Astrophys.J .Suppl. 180 
(2009) 330-376, [arXiv: 0803. 0547]. 

[35] C. Burgess, R. Easther, A. Mazumdar, D. F. Mota, and T. Multamaki, Multiple inflation, 
cosmic string networks and the string landscape, JHEP 0505 (2005) 067, [hep-th/0501125]. 

[36] K.-Y. Choi, S. A. Kim, and B. Kyae, Primordial curvature perturbation during and at the end 
of multi-field inflation, arXiv: 1202 .0089. 20 pages. 

[37] C. T. Byrnes, K.-Y. Choi, and L. M. Hall, Large non-Gaussianity from two- component hybrid 
inflation, JGAP 0Q02 (2009) 017, [arXiv: 0812. 0807]. 

[38] Q.-G. Huang, A Geometric description of the non-Gaussianity generated at the end of 
multi-field inflation, JC^P 0906 (2009) 035, [arXiv: 0904. 2649]. 19 pages, 4 figures. 

[39] R. AUahverdi, K. Enqvist, J. Garcia-Bellido, and A. Mazumdar, Gauge invariant MSSM 
inflaton, Phys. Rev. Lett. 97 (2006) 191304, [hep-ph/0605035]. 

[40] R. AUahverdi, A. Kusenko, and A. Mazumdar, A-term inflation and the smallness of neutrino 
masses, JCAP 0707 (2007) 018, [hep-ph/0608138]. 

[41] R. AUahverdi, K. Enqvist, J. Garcia-Bellido, A. Jokinen, and A. Mazumdar, MSSM flat 
direction inflation: Slow roll, stability, fine tunning and reheating, JGAP 0706 (2007) 019, 



20 



[hep-ph/0610134]. 

[42] S. Clesse, C. Ringeval, and J. Rocher, Fractal initial conditions and natural parameter values in 
hybrid inflation, Phys.Rev. D80 (2009) 123534, [arXiv: 0909. 0402]. 

[43] F. Duplessis, Y. Wang, and R. Brandenberger, Multi-Stream Inflation in a Landscape, JCAP 
1204 (2012) 012, [arXiv: 1201.0029]. 

[44] S. Dimopoulos, S. Kachru, J. McGreevy, and J. G. Wacker, N-flation, JCAP 0808 (2008) 003, 
[liep-th/0507205]. 

[45] R. Easther and L. McAllister, Random matrices and the spectrum of N-flation, JCAP 0605 
(2006) 018, [hep-th/0512102]. 

[46] S. Yokoyama, T. Suyama, and T. Tanaka, Primordial Non- Gaussianity in Multi-Scalar 
Slow- Roll Inflation, JCAP 0707 (2007) 013, [arXiv: 0705. 3178]. 

[47] S. Yokoyama, T. Suyama, and T. Tanaka, Primordial Non- Gaussianity in Multi-Scalar 
Inflation, Phys.Rev. D77 (2008) 083511, [arXiv: 0711 .2920]. 



21 



